Ballistic versus diffusive magnetoresistance of a magnetic point contact 
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The quasiclassical theory of a nanosize point contacts (PC) 
between two ferromagnets is developed. The maximum avail- 
able magnetoresistance values in PC are calculated for ballis- 
tic versus diffusive transport through the area of a contact. 
In the ballistic regime the magnetoresistance in excess of few 
hundreds percents is obtained for the iron-group ferromag- 
nets. The necessary conditions for realization of so large mag- 
netoresistance in PC, and the experimental results by Garcia 
et al are discussed. 

PACS numbers: 74.80. Dm, 74.50.+r, 74.62.-c 



In recent experiments on study of Ni-Ni and Co-Co 
point contacts (PC), a surprisingly high negative m. 
netoresistance exceeding 200% has been discovered 
The set up of the experiment was typical for observation 
of giant magnetoresistance (GMR), the effect observed 
earlier in hybrid systems, involving ferromagnetic and 
normal multilayer metals .EHj However, for the multilayer 
structures the typical change of the resistance reached 
10% — 30%, which is considerably lower than the corre- 
sponding values of Refs. |l|,||. So, one can come easily 
to the conclusion that the main contribution to the MR 
comes from the region of the PC. 

A negative magnetoresistance can be due to scattering 
on domain walls (DW) ancLthis effect has been consid- 
ered in a number of worksOu! giving typical values of MR 
in a range of few percents. Such considerably low values 
of the MR were obtained assuming that realistic widths 
of the DW were large, which resulted in a low scatter- 
ing amplitudes. Considering sharrtDW in the ballistic 
regime one comes to values ~ 70%Ja 

The fact that a sharp DW may give large MR was used 
in Ref. || to explain the anomalously large. values of MR 
in the experiments on the point contacts.Era However, the 
theory in Ref. 2 is the perturbation theory, it can not be 
applied to the explanation of 300% effect. The dimin- 
ishing of the width of DW, when decreasing nfche size of 
the constriction was demonstrated by Bruno The DW 
width becomes comparable with PC length, and magne- 
tization rotates almost abruptly inside the constriction. 
This conclusion holds until the diameter of PC is smaller 
than its actual length. With further increase of constric- 
tion size (diameter) the wall will bend outside of PC, and 
simple energy considerations ahpw that the DW width 
will be of the order of PC sizc.E3 

The regime of conductance quantization in magnetic 



PC has been considered by Imamura et al .0 They ob- 
tained that, if the spin of conduction electron cannot 
rotate in DW pinned to the constriction, then magne- 
toresistance acquires oscillations as a function of PC size 
with amplitude exceeding 1000%. 

In this paper we develop a quasiclassical theory of 
electric transport through magnetic PC taking into ac- 
count scattering by impurities, thus covering the ballistic 
(I > a) and diffusive (I < a) regimes (I is the mean free 
path and a is the radius of the contact). The typical PC 
size, which is beyond the quantization regime, 2a > 8A, 
may be well described within the quasiclassical (QC) ap- 
proximation (2a 3> Aj? = 2-7T /pf ~ 6 A, Xp and pf are 
the Fermi wavelength and momentum). 

We believe that extremely large magnetoresistance can 
be obtained, if the strong reflection of spin-polarized cur- 
rent carriers on the PC area is achieved at antiparallel 
(AP) alignment of magnetizations in contacting ferro- 
magnets. This is realized, if there is mismatch in the spin- 
subbands Fermi- momenta of contacting magnets. For AP 
alignment pfit = PF2[, and pfii — PF2]- Let us assume 
that pfit 3* Pfii- Then subband with the smaller value 
of the Fermi momenta, which is minority subband, can 
not accept momenta transferred from the opposite side 
of the PC, which is majority subband with the same spin 
projection. As a result, only a narrow incidence angles 
cone around the normal to the interface is responsible for 
the charge transport across the PC. Electrons with more 
inclined trajectories are completely reflected. Thus, the 
partial transmission at the steep incidence, and the to- 
tal reflection at slanting incidence provide high boundary 
resistance of PC. 

The necessary condition for realization of the above 
scenario is the conservation of electron spin orientation 
when crossing the domain wall. The orientation con- 
serves, if the DW width d w is shorter than the length 
d s , at which the electron spin quantization axis adjusts 
the varying direction of local exchange field. For bal- 
listic transmission through PC d s — vfT\, where T\ is 
the longitudinal relaxation time of conduction electron 
magnetization - the Overhauser time.t 2 ! At this condition 
the transmission process looks like transmission through 
abrupt DW, and the description of the electron transport 
through PC with boundary conditions at PC interface is 
valid. 

The PC model we consider is the circular hole of the 
radius a made in a membrane, which divides the space 
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on two half-spaces, occupied by single-domain ferromag- 
netic metals. The membrane is impenetrable for the 
quasiparticles carrying a current, however the thickness 
of the membrane in the model is assumed to be vanish- 
ing. The z— axis of coordinate system is chosen perpen- 
dicular to the membrane plane. The electron motion on 
both sides of the contact can be described by the equa- 
tions fowjuasiclassical (QC) Green functions derived by 
ZaitsevB They are in fact the Boltzmann equations in 
the r— approximation: 



dg a . dg s 1 



= 0, 



dg s 
dz 



dga 



^ = 0. 



(1) 



g s and g a are symmetric and antisymmetric with re- 
spect to z— projection of quasiparticle momentum QC 
GF (Green functions integrated over the energy vari- 
able), v is the vector of the Fermi velocity, v z — vp cos 6, 
v f\ = v 'f ~~ v z i angle 9 is measured from the z— axis, vf 
is the modulus of v, the bar over ~g c means the averag- 
ing over the solid angle. We assume that the spin-mixing 
process is weak, therefore we consider spin channels as in- 
dependent and omit the spin-channel indices in (Q) and 
expressions below. 

The boundary conditions to equation ([!]) for the spec- 
ular scattering ( Pf± a shi0i = PF2aSin#2 = at the 
interface z = are:E3 

/ n N fn , J 3a(0), p\\ <PF1,PF2 

1^0, mm(p F i,p F2 ) < p\\ 



2Rg a (0) = -D(g s2 -g sl ), 



(2) 



where subscript 1 or 2 labels left- or right-hand side of 
the contact, respectively, p F i is the Fermi momentum of 
i-th side, p|| is the projection of the Fermi momentum 
vector on the PC plane. D and R= 1 — D are the exact 
quantum mechanical transmission and reflection coeffi- 
cients that can be considered either as phenomenological 
parameters or calculated for models of interest. The sec- 
ond line in the first boundary condition in ( |^) explicitly 
quantifies the total reflection for inclined trajectories, de- 
scribed qualitatively above. 

The density of a current through the contact may be 
written as 
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dfte cos 6 <g a (z, ~p*, t) . (3) 



The total current through the area of the contact is 

oo 

I z (z -> 0,i) = a / dkJi{ka)f(0,k,t). (4) 
o 

In the above equations pFmin = min(pFi,PF2), Ji(x) is 
the Bessel functions, j z (0,k,t) is the Fourier-transform 



of current density, Eq. (||), over the in-plane coordinate 
p . The cylindrical symmetry of the problem has been 
used upon derivation of Eq. (^) . 

We search a solution for g s in the form (Ub = ft = 1): 



g s (e) = tanh — + / s (e) 



(5) 



where the first term is the equilibrium value of g s in the 
leads far away of PC. Substitution of @ into (Q) and 
Fourier transformation over the variable p leads to equa- 
tions, the exact solution of which reads 



f s (z) = 9a (z)sgn(z) + f- / d£e-^-*l/ s (£,fc), (6) 



where 



1 - ikl 



(7) 



I = rvp is the mean free path, l z = IcosO, ijj = I 2 — l\. 
Integrating Eq. (||) over the solid angle we obtain 



f s (z>0)=g a + dZK(Z~z)f s (Z,k), 



(8) 



where the kernel K{rf) is (x = cos ( 
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If the mean free path I is short (I -C a), the second term 
in Eqs. (^) and (@) dominates and the integrand of Eq. 
(Q) is the product of rapidly decreasing on the distance 
I kernel K{rj) and slowly decreasing function f s . That 
is why we may take out f c (k,£) from the integral (|§|) at 
the point £ = z. Within this approximation we obtain 



f.(z,k)=g a (z,k)(l-X(k)Y 



(10) 



where 



X(k) = fd£KU-z) = iarctanJW. (11) 
J kl 



Substituting Eq. ( |10[ ) into Eq. (|6|) , and using the bound- 
ary conditions (0) we obtain the equation for the anti- 
symmetric combination g a : 



9a(0,k) 



1 r,f , e , e-eV\ 

—D tanh tanh 

2 \ 2T 2T J 



Ik 



D 



1 



1 - Ai 2l zl 
D 1 

1 — A2 2l z 2 



d^g al (0 
dHe~^y- a2 (0, (12) 



2 



where 
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= I pdp J <' k// </ 7 ; 
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J\ (ka) 



(13) 



V is the bias voltage. 

To find g a (0,k) we average Eq. ( |l2] ) over the solid 
angle, exploit the continuity of g a at the interface, Eq. 
(^), and, again, use the fact that in the limit I <C a the 
kernel in the integral over x in the second and third terms 
of averaged Eq. ( |l2] ) is a function, rapidly decreasing 
at distance /. Of course, in the ballistic regime (I > 
a ) this approximation is no longer valid, but, in this 
regime the first (exact) term in Eq. ([l2]) dominates the 
approximate terms with integrals. So, in the ballistic 
limit the approximation done does not bring a big error 
either. Although the approximation may not be valid 
in the intermediate regime, the suggested scheme can be 
used as an interpolation. 

Now we find easily <L(0, k), make consecutive substi- 
tutions into ( |l2| ) , (||) , ([|) , and, finally, obtain the general 
expression for the current through PC: 



e PF min a V fdk 
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where 

F{k,6) = 1 - 



-Jl(ka) (DF(k, 9) cos 0) , (14) 
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2 (1 - Ai) xtUi 2(l-A 2 )x 2 Z* 2 
D 



i+ 2(1-Ax) + 2(1-A 2 ) 



(15) 
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^/l + k 2 l 2 {l-x 2 Y 



(16) 



< . . . > means averaging over the solid angle. Eqs. (0) 
and ( |l5| ) are the basic analytical result of the paper, 
which expresses the current in terms of parameters D, 
I, a, pf characterizing the system. 

Now we calculate the magnetoresistance of PC between 
two identical ferromagnets. It can be expressed via the 
conductances a — I/V as follows 



MR = 



R 



AP 



R 1 



r AP 



R p 



rAP 



(17) 



where R p (<r p ) stands for the resistance (conductance) 
at parallel alignment of magnetizations of contacting fer- 
romagnets, and R AP (cr AP ) is for the antiparallel align- 
ment of magnetizations. For the parallel alignment the 
net current is the sum of currents for both (independent) 
spin channels, D — 1, Ai — Ai. Labelling the quantities 
by arrow- up/down notations we write down 



= (T tT +<7 U 



' 2 (ph+P 2 Fl) (™ 2 ) Jdk 
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J t jKka) 



ph 



+ (T^l) } ■ (18) 



Pp 1 +Pp l (l + Jl + k 2 l 2 



The prefactor in Eq. ( Jig ) is nothing but the sum of 
SharvinE-a conductances for the spin channels. For the 
AP alignment of magnetizations the conductance is 
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2 (1 - AT) (1 - Ai) + Aj^l - Ai) + A| x (l - AT) J 

where (D(x))^ stands for the transmission coefficient 
of the interface at AP alignment. For the mechanism 
of magnetoresistance discussed above, (D(x)) + , can be 
found from the solution of Schrodinger equation fop-jthe 
particle moving in the step-like potential landscape^ 



(19) 
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with v z2 = v zl for the antiparallel alignment. The trans- 



mission coefficient (|20|) gives maximum available mag- 
netoresistance values for a particular parameters choice. 
Neglecting the difference of the effective masses in the 
spin-subbands we may write down 



(D(x)) 



Ax\fb 2 



U 



y/b 2 + x 2 ) 



where 



(21) 
(22) 

For the purely ballistic transport (a/Zj — > 0, where 1^ 
(Zj) is the majority (minority) electrons mean free path) 
all integrals in Eqs. ©, © 

are evaluated analytically, 

and magnetoresistance reads 



S 2 



PFl 
PF\ 



< 1. 



MR 



{1 - 5) {56 3 + 155 2 + 95 + 3} 
8S 3 {5 + 2) ■ 



(23) 



If 5 = 1, then MR = 0, i.e. the magnetoresistance 
vanishes. For the set of S values we obtain from (E2 
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S = 0.5, MR = 238%; S = 0.4, MR = 455%; S = 0.33, 
MR = 780%; <5 = 0.3, MR = 1012%. 

In general case the angular integrals in (|l6|) and jl9| ) 
can be still evaluated analytically, whereas the integra- 
tions over k can be done only numerically. The results for 
the magnetoresistance ( |l7|) as function of the contact ra- 
dius are shown on Fig. 1. The curves show the maximum 
available MR, that could be realized in PC with physical 
parameters displayed on the figure. MR exponentially 
drops when the size of the contact approaches the mean 
free path of a material. Then it shows a smooth crossover 
from ballistic to diffusive regimes of conduction. 

Let us discuss the experimental data on magnetore- 
sistance of magnetic PC by Garcia et al. Ni-Ni PC 
showed maximal MM ~ 280%,H Co-Co PCUjhowed max- 
imal MR ~ 230%E In the recent paperB they quote 
maximal MR ~ 33% for Fe-Fe PC. To obtain the MR 
values 280% (Ni) and 230% (Co) we have to use the 
values (5(Ni) ~ 0.47 and <5(Co) ~ 0.5. These numbers 
are in the range of the values, obtained experimentally 
frortuthe single photon threshold photoemission, 8{G6) ~ 
0.4,t2l and from ferromagnet/supercond-uctor point cop-, 
tact spectroscopy: SM) ~ 0.59 - 0.65X3 5(Ni) ~ 0.71 ;ta 
<J(Co) ~ 0.62 - 0.65,0 S(Co) ~ 0.68.E3 

If we use the experimental data of Ref. [Tt] for iron, 
<5(Fe) ~ 0.59 — 0.65, then in our theory we obtain 
M.R(Fe) = (100- 140)% r which is larger than the experi- 
mentally measured 33%Jl3 The justification of our model 
suggests that observed MR does not solely confined to a 
value of polarization S. We believe that the basic condi- 
tion for observation of upper MR limit, d w <C d s , is not 
fulfilled in the Fe-Fe PC experiment £3 T\ is proportional 
to the squared magnetic moment and the integral of ex- 
change between conduction electrons and localized mo- 
ments, and proportional to conduction electrons density 
of states at Fermi level. All these physical parameters 
for iron are larger than for cobalt, and especially than 
for nickel. Therefore we expect, that Ti(Fe) at least one 
order of magnitude shorter than T\ for Co and Ni. When 
d s (Fe) ~ d w (Fe), the electron spin almost tracks the lo- 
cal exchange field in the domain wall. As a result the 
reflection of the electrons from DW decreases, and the 
observed MR does not reach its maximal value. 

Let us discuss now the magnetoresistance in the diffu- 
sive regime of transport, when the radius of nanohole 
is much larger than the mean free path (a ^> 
The giant MR values can be obtained, if the condition 
of validity of our model, d w <C d s , will be realized in an 
experiment. In the opposite limit, d w > d s , when PC 
size is so large that DW becomes smooth and wide, the 
electron spin will track the local exchange fieldJn the do- 
main wall, and MR will level off at Levy-Zhangu impurity 
scattering enhancement mechanism, which can give 2- 
11% magnetoresistance. The requirement of abrupt DW 
with constant width, irrespective of the PC size, can be 
technologically controlled, if very thin (2-4 monolayers of 
the thickness ~ Ap) nonmagnetic interlayer is deposited 
on the PC plane before depositing the second electrode. 



Then, just like in CPP transport in multilayers ,trH the 
contacting domains will be exchange decoupled, so the 
magnetization will acquire sudden reversal within the 
spacer thickness ~ Xf- In this case our analysis is valid 
for an arbitrary size of PC. 
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Figure captions 

Fig.l. The dependence of magnetoresistance on the 
PC radius. 
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